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ABSTRACT:  ’/-Glutamyl  transpeptidase  (GGT)  is  a  two- 
substrate  enzyme  that  plays  a  central  role  in  glutathione 
metabolism  and  is  a  potential  target  for  drug  design.  GGT 
catalyzes  the  cleavage  of  /-glutamyl  donor  substrates  and  the 
transfer  of  the  /-glutamyl  moiety  to  an  amine  of  an  acceptor 
substrate  or  water.  Although  structures  of  bacterial  GGT  have 
revealed  details  of  the  protein— ligand  interactions  at  the  donor 
site,  the  acceptor  substrate  site  is  relatively  undefined.  The 
recent  identification  of  a  species-specific  acceptor  site  inhibitor, 

OU749,  suggests  that  these  inhibitors  may  be  less  toxic  than 
glutamine  analogues.  Here  we  investigated  the  donor  and 
acceptor  substrate  preferences  of  Bacillus  anthracis  GGT 
(CapD)  and  applied  computational  approaches  in  combina¬ 
tion  with  kinetics  to  probe  the  structural  basis  of  the  enzyme’s  substrate  and  inhibitor  binding  specificities  and  compare  them 
with  human  GGT.  Site-directed  mutagenesis  studies  showed  that  the  R432A  and  R520S  variants  exhibited  6-  and  95-fold 
decreases  in  hydrolase  activity,  respectively,  and  that  their  activity  was  not  stimulated  by  the  addition  of  the  L-Cys  acceptor 
substrate,  suggesting  an  additional  role  in  acceptor  binding  and/or  catalysis  of  transpeptidation.  Rat  GGT  (and  presumably 
HuGGT)  has  strict  stereospecificity  for  L-amino  acid  acceptor  substrates,  while  CapD  can  utilize  both  l-  and  D-acceptor 
substrates  comparably.  Modeling  and  kinetic  analysis  suggest  that  R520  and  R432  allow  two  alternate  acceptor  substrate  binding 
modes  for  l-  and  D-acceptors.  R432  is  conserved  in  Francisella  tularensis,  Yersinia  pestis,  Burkholderia  mallei,  Helicobacter  pylori  and 
Escherichia  coli,  but  not  in  human  GGT.  Docking  and  MD  simulations  point  toward  key  residues  that  contribute  to  inhibitor  and 
acceptor  substrate  binding,  providing  a  guide  to  designing  novel  and  specific  GGT  inhibitors. 


/-Glutamyl  transpeptidase  (GGT,  EC  23.2.2)  is  an  enzyme 
that  plays  a  central  role  in  glutathione  metabolism  and  cysteine 
salvage.1  3  It  has  been  found  in  a  wide  variety  of  organisms, 
including  mammals,  plants,  and  bacteria.  The  latter  includes 
such  pathogens  as  Helicobacter  pylori,  Escherichia  coli, 
Burkholderia  mallei,  Staphylococcus  aureus,  Yersinia  pestis,  Bacillus 
anthracis,  and  Francisella  tularensis.1’2  In  humans,  GGT 
(HuGGT)  is  implicated  in  a  number  of  physiological  processes 
and  disorders,  such  as  tumor  progression,  drug  resistance,  and 
neurodegenerative  diseases.4  The  bacterial  GGT  enzyme  in 
H.  pylori  (HpGGT)  is  an  apoptosis-inducing  protein  and  is 
essential  for  the  establishment  of  infection.7’8  The  B.  anthracis 
GGT  (CapD)  is  a  virulence  factor  required  in  anchoring  the 
capsule  to  the  bacterial  cell  wall  as  well  as  capsule  remodeling.9 
CapD  is  an  unusual  GGT  in  that  it  has  been  shown  to  be 
unable  to  hydrolyze  glutathione  (GSH);3  its  endogenous 
substrate  is  likely  the  poly-y-D-glutamic  acid  (PDGA)  capsule 
surrounding  the  bacterium.  The  capsule  allows  B.  anthracis  to 
evade  the  immune  system  by  preventing  phagocytic  killing  by 
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immune  cells  and  by  acting  as  a  physical  barrier  to  surface 
antigens.  The  capsule  itself  is  weakly  immunogenic.  CapD  is 
not  a  typical  antimicrobial  target  in  that  it  is  not  a  bactericidal 
target.  However,  inhibition  of  CapD  results  in  the  loss  of  the 
capsule  from  the  cell  wall,  resulting  in  unencapsulated  bacteria 
that  are  now  susceptible  to  phagocytosis  by  neutrophils  and  can 
be  cleared  by  the  immune  system.  Moreover,  CapD  is  a 
membrane-associated  enzyme,3  and  inhibitors  of  the  enzyme 
could  be  used  to  bypass  some  mechanisms  of  antimicrobial 
resistance  (e.g.,  efflux  pumps).  GGTs  of  Gram-negative  bacteria 
are  also  potential  antimicrobial  targets  as  host  glutathione  can 
function  as  a  source  of  cysteine.  Mutant  strains  of  F.  tularensis 
that  lack  functional  GGT  have  been  shown  to  have  impaired 
intracellular  replication  and  are  strongly  attenuated  for 
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virulence  in  mice.11  Therefore,  GGT  has  been  identified  as  a 
valid  target  for  both  anticancer  and  antimicrobial  therapeutics. 

GGTs  catalyze  the  cleavage  of  /-glutamyl  donor  substrates 
(e.g.,  glutathione  and  its  S-conjugates),  and  the  transfer  of 
the  /-glutamyl  group  to  an  acceptor  substrate  amino  acid, 
N-terminal  residue  of  a  polypeptide  (transpeptidation),  or  water 
(hydrolysis)  following  a  ping-pong  kinetic  mechanism.1”’13 
GGT  is  expressed  as  an  ~60  kDa  precursor  and  is  activated  by 
an  autocatalytic  internal  cleavage  (e.g.,  S  jTTHFV  in  CapD  and 
GjTAHLS  in  HuGGT)  to  generate  a  heterodimer  consisting 
of  ~40  and  ~20  kDa  subunits.  The  N-terminal  threonine 
(underlined)  hydroxyl  of  the  small  subunit  acts  as  the 
nucleophile  in  both  the  autoprocessing  and  enzyme-catalyzed 
reactions.1'" 

A  close  homologue  of  HuGGT,  namely,  rat  GGT  that  is  79% 
identical  to  human  GGT,  can  utilize  a  number  of  acceptor 
substrates  in  a  pH-dependent  fashion.  Typical  transpeptidase 
acceptor  substrates  have  a  terminal  amino  group.  The  simplest 
acceptor  substrates  include  single  amino  acids  and  glycine 
dipeptides  (Gly-Ala  or  Z-Gly,  where  Z  is  W,  F,  Y,  H,  A,  Q,  C, 
M,  S,  L,  K,  A,  or  G).13  The  product  of  the  GGT-catalyzed 
hydrolysis  of  glutathione,  Cys-Gly,  mimics  both  the  donor 
substrate  and  dipeptide  acceptor  substrates.  The  acceptor 
substrate  can  therefore  exploit  features  of  the  binding  site  that 
stabilize  binding  of  glutathione  on  the  acceptor  side.  Cook  et  al. 
recognized  that  the  pXa  of  the  a-amino  group  of  glycyl 
dipeptides  is  significantly  lower  than  those  of  free  amino  acids 
and  may  in  part  explain  why  glycyl  dipeptides  tend  to  be  better 
acceptor  substrates  for  rat  GGT  than  free  amino  acids.  The 
transpeptidation  reaction  has  a  bell-shaped  pH— rate  profile  for 
both  rat  GGT  and  HuGGT,14'16-18  but  only  an  ascending  limb 
in  E.  coli  GGT.11  The  hydrolysis  reaction  yields  a  relatively  flat 
pH— rate  profile.  While  free  amino  acids  are  typically  found  in 
their  zwitterionic  forms,  the  reactive  form  of  the  acceptor 
substrate  is  believed  to  contain  a  neutral  terminal  amine  as  the 
pJCa  of  the  ascending  limb  varies  with  the  acceptor.15  On  the 
acceptor  substrate  side,  a  positively  charged  residue  with  a  pK, 
of  8.46  in  rat  GGT  (a  similar  pH— rate  profile  was  reported  for 
HuGGT18)  was  proposed  to  be  important  for  binding  the 
acceptor  and/or  catalysis;15  however,  a  specific  residue  has  not 
been  shown  to  account  for  the  descending  limb  of  the  pH— rate 
profile.  Herein,  we  identify  an  important  mobile  electrostatic 
catalyst,  Arg-520  in  CapD  (Lys-562  in  HuGGT),  which  may 
interact  with  the  acceptor  substrate  and  contribute  to  catalysis. 

X-ray  structures  of  GGTs  from  several  bacterial  species  have 
revealed  atomic  details  of  protein— donor  substrate  interac- 
tions  "  and  the  acyl  intermediate.  Despite  a  high  degree  of 
conservation  of  structure  and  function  between  human  and 
bacterial  GGT  homologues,  significant  differences  in  acceptor 
substrate  and  inhibitor  preferences  are  evident.  For  example,  rat 
and  E.  coli  GGT  have  strict  stereospecificity  for  L-amino  acid 
acceptor  substrates,26’27  while  B.  anthracis  CapD  does  not  (this 
paper).  HuGGT  also  catalyzes  transpeptidation  180-fold  faster 
than  hydrolysis,28'21  while  HpGGT  shows  a  modest  2-fold30 
and  CapD  a  12-fold  difference  (this  paper)  between  the  rates. 
Residues,  which  account  for  the  faster  rate  transpeptidation 
versus  hydrolysis,  and  acceptor  substrate  stereospecificity  have 
not  been  previously  identified.  We  hypothesized  that  differ¬ 
ences  may  arise  from  nonconserved  residues  in  the  relatively 
undefined  acceptor  substrate  binding  site. 

Within  the  putative  acceptor  substrate  site,  two  alternate  side 
chain  conformers  were  observed  for  Arg-520  in  the  X-ray 
crystal  structure  of  CapD  [Protein  Data  Bank  (PDB)  entry 


3G9K25].  To  determine  if  the  flexibility  of  this  and  other 
residues  may  also  account  for  differences  in  substrate  binding 
and/or  product  release,  we  examined  molecular  dynamics 
simulations  of  bacterial  and  human  GGTs. 

Human  GGT  is  a  complementary  drug  target  for  various 
glutathione-sensitive  cancer  therapies.  ’1-  ’4  Although  a  number 
of  potent  inhibitors  of  HuGGT  have  been  discovered  in  the 
past,  most  modify  the  enzyme  irreversibly  and  have  associated 
toxicities.  Richter  et  al.  also  identified  an  inhibitor  for  CapD, 
Capsidin,9  which  irreversibly  acylates  the  active  site  threonine. 
While  effective  in  inhibiting  the  enzyme  in  vitro,  the  acetyl 
groups  of  Capsidin  were  readily  hydrolyzed  by  serum 
carboxylesterases.  The  human  GGT  inhibitors,  which  bind  in 
the  donor  substrate  site  and  mimic  glutamate,  such  as  acivicin, 
inhibit  CapD  poorly  (ICS0  >  20  mM)  (data  not  shown)  and 
cannot  be  repurposed  as  antimicrobials.  Recently,  a  novel  class 
of  HuGGT  inhibitor  OU749  was  identified.34  Unlike  donor 
substrate  analogue  inhibitors,  OU749  is  believed  to  bind  the 
acceptor  site  and  is  an  uncompetitive  inhibitor  with  respect  to 
the  /-glutamyl  donor  substrate.  Consequently,  OU749  is 
significantly  less  toxic  than  acivicin  and  exhibits  strong 
species-specific  GGT  inhibition.34  The  discovery  of  OU749 
suggested  a  promising  direction  toward  the  design  of  nontoxic, 
species-selective,  acceptor  site  inhibitors  of  GGT,  which  may 
extend  to  bacterial  GGT. 

In  this  study,  we  applied  computational  and  kinetic 
approaches  to  investigate  both  the  donor  and  acceptor 
substrate  specificities  of  B.  anthracis  CapD  and  compare  them 
to  those  of  HuGGT  to  identify  any  unique  features  of  the 
bacterial  enzyme’s  active  site  that  could  be  exploited  in  follow- 
on  in  silico  drug  screening.  In  addition  to  natural  acceptor  and 
donor  substrates  and  their  derivatives,  we  examined  the 
uncompetitive  acceptor  site  inhibitor  OU749  as  a  molecular 
probe  to  further  explore  the  structural  basis  of  inhibitor 
binding.  One  surprising  commonality  between  rat  GGT  (and 
presumably  HuGGT,  which  is  79%  identical)  and  CapD  is  that 
both  GGTs  share  the  same  acceptor  substrate  (L-Cys) 
preference  (this  paper)  but  differ  in  their  binding  specificity 
for  OU749.  To  begin  to  understand  species-specific  inhibition 
within  the  acceptor  substrate  site,  we  determined  the  optimal 
binding  conformation  of  selected  substrates  and  inhibitors 
using  a  multistage  docking  and  refinement  strategy.  Molecular 
dynamic  simulations,  site-directed  mutagenesis,  and  kinetic 
analysis  were  used  to  characterize  key  residues  at  the  putative 
acceptor  substrate  site  that  may  account  for  species-specific 
acceptor  substrate  and  acceptor  site  inhibitor  binding 
specificities. 

■  EXPERIMENTAL  PROCEDURES 

B.  anthracis  CapD  Cloning,  Protein  Expression,  and 
Purification.  The  CapD  expression  construct  expressed  a 
maltose-binding  protein  (MBP)  fusion  protein  with  a  Factor  Xa 
cleavage  site  prior  to  the  CapD  sequence.  The  first  27  residues 
were  predicted  to  be  the  signal  peptide  sequence  and  were 
excluded.35  The  insert  was  placed  between  the  XmnI  and 
Hindlll  cleavage  sites  of  the  pMAL-c2  plasmid  (New  England 
Biolabs,  Ipswich,  MA)  using  the  forward  primer,  5'CAGC- 
GAAGGATTTCTTTCAATAAAATAAAAGACAGTGTTA3', 
and  the  reverse  primer,  5'CAGCAAGCTTCTATTTATTT- 
GATTTCCAAGTTCCATT3'.  The  protein  was  expressed  in 
E.  coli  BL-2l(DE3)  cells  at  17  °C  for  20  h.  Cells  were  lysed  in 
IX  phosphate-buffered  saline  (PBS)  (pH  7.4),  5%  glycerol, 
0.3  mg/ mL  lysozyme,  and  0.05  mg/mL  DNAase  and  sonicated. 
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The  lysate  was  centrifuged  for  30  min  at  20000g,  and  the 
supernatant  was  loaded  onto  an  amylose-resin  column 
equilibrated  with  IX  PBS  (pH  7.4).  The  column  was  eluted 
with  IX  PBS  (pH  7.4)  and  10  mM  maltose.  The  protein 
was  dialyzed  into  0.5X  PBS  (pH  7.4)  and  loaded  onto  a 
CTSepharose  column.  The  protein  was  eluted  with  a  gradient  of 
sodium  chloride  (0  to  1.0  M).  Fractions  were  combined  and 
dialyzed  against  0.5X  PBS  (pH  7.4)  and  loaded  onto  an  SP- 
Sepharose  column  equilibrated  with  the  same  buffer.  The  pro¬ 
tein  was  eluted  with  a  gradient  of  sodium  chloride  (0  to  1.0  M). 
The  protein  was  concentrated,  and  the  buffer  was  exchanged 
into  IX  PBS  (pH  7.4)  and  50%  glycerol  and  stored  at  —20  °C. 
For  assays,  the  glycerol  was  removed  using  a  PD-10  column 
equilibrated  with  assay  buffer. 

Expression  of  B.  anthracis  CapD  Protein  and 
Purification  of  Site-Directed  Mutants.  Three  site-directed 
mutants,  R520S,  R432A,  and  N373A,  were  constructed  by 
Geneart  Inc.  (Regensburg,  Germany)  using  the  pMAL  plasmid. 
All  were  purified  using  the  same  procedure  as  for  wild-type 
(WT)  CapD-MBP  described  above. 

B.  anthracis  CapD  Assays.  Capsule  degradation  was 
observed  by  gel  electrophoresis  using  a  0.5  to  1.0%  agarose 
gel  in  IX  TAE  buffer  (40  mM  Tris-acetate  and  1  mM  EDTA). 
The  capsule  was  purified  from  the  Ames  strain  of  B.  anthracis  as 
described  previously.36  The  gel  was  stained  with  Loeffler’s 
methylene  blue  solution  (0.3  g  of  methylene  blue,  30  mL  of 
ethanol,  and  100  mL  of  water).  Reaction  mixtures  contained 
100  mM  Tris  (pH  8.0),  capsule,  CapD-MBP  (25  fig/mL),  and 
acceptor  (0.5  /^g/ mL)  and  were  incubated  for  30  min  at  37  °C. 
The  capsular  material  was  used  at  a  final  concentration  of  0.8— 
1.0  fig/fih  for  assays  and  was  held  at  a  constant  concentration 
in  an  experiment. 

Steady  state  kinetic  parameters  were  measured  using  a 
synthesized  FRET  donor  substrate,  Abz-[y-D-Glu]5-Dnp, 
described  by  Richter  et  al.9  The  FRET  donor  substrate  is  a 
5-mer  of  poly-D-glutamate  with  an  N-terminal  2-aminobenzoyl 
(Abz)  fluorophore  and  C-terminal  2,4-dinitrophenyl  (DNP) 
quencher.  The  donor  substrate  was  synthesized  by  BioPeptide 
Co.  Inc.  (San  Diego,  CA).  Initial  velocities  were  measured  by 
varying  the  concentration  of  the  FRET  donor  substrate  and 
acceptor  substrate.  Similar  kinetic  analyses  of  GGT  have  been 
conducted  by  others.26,37  Initial  velocities  were  globally  fit  to 
the  velocity  equation  for  a  ping-pong  bi-bi  kinetic  mechanism: 

ymJS][A] 

Kma[S]  +Kms[A]  +  [S][A]  (1) 

where  [S]  is  the  concentration  of  the  FRET  donor  substrate, 
[A]  is  the  concentration  of  the  acceptor  substrate,  K ma  is  the 
Michaelis  constant  for  the  acceptor  substrate,  and  Kms  is  the 
Michaelis  constant  for  the  donor  substrate.  The  Kms  of  the 
donor  substrate  is  an  apparent  Km  as  the  acceptor  can 
competitively  inhibit  the  enzyme.13,26  Rates  were  measured  in 
25  mM  HEPES  (pH  7.4)  and  0.1%  Tween  20  at  22  ±  3  °C 
using  an  excitation  wavelength  of  315  nm  and  an  emission 
wavelength  of  415  nm.  Data  were  corrected  for  inner  filter 
effects  using  the  same  substrate  lacking  the  quencher.  Data 
were  fit  using  Grafit  version  6.0.1  (Erithacus  Software  Ltd.).  All 
chemicals  were  purchased  from  Sigma  (St.  Louis,  MO). 

Expression  and  Purification  of  HuGGT.  The  HuGGT 
expression  vector  was  a  kind  gift  from  M.  Hanigan  (University 
of  Oklahoma,  Oklahoma  City,  OK).  COS-1  cells  were  plated  in 
three  75  cm2  tissue  culture  flasks  in  10  mL  of  Dulbecco’s 


modified  minimum  essential  medium  (DMEM)  (Cellgro) 
containing  10%  fetal  bovine  serum  (Cellgro)  at  37  °C  in  the 
presence  of  5%  C02.  The  following  day,  the  flasks,  at 
approximately  80%  confluence,  were  transfected  with  5.5 
of  the  pcDNA3.l(+)  HuGGT  plasmid  and  20  ftL  of 
Lipofectamine  2000  (lipofectamine)  (Life  Technologies), 
according  to  the  manufacturer’s  recommendations.  Following 
a  48  h  incubation,  the  cells  were  washed  twice  with  ice-cold  IX 
PBS  and  then  harvested  by  being  scraped.  The  cells  were 
centrifuged  at  1000  rpm  and  4  °C  for  15  min,  and  the  pellet 
was  flash-frozen  in  liquid  nitrogen.  Cells  were  lysed  in  1.5  mL 
of  lysis  buffer  [50  mM  HEPES  (pH  7.4),  150  mM  NaCl,  0.5% 
CHAPS,  and  1  EDTA-free  protease  tablet  (Roche)  per  2.5  mL 
of  lysis  buffer]  for  30  min  with  rocking.  Protein  was  dialyzed 
into  500  mL  of  50  mM  HEPES  (pH  8.0),  150  mM  NaCl,  with 
1.0  mL  of  uncharged  Chelating  Sepharose  (GE  Healthcare)  to 
remove  any  metal  ions  and  the  detergent.  Protein  was  loaded 
onto  a  Chelating  Sepharose  column  (GE  Healthcare)  charged 
with  nickel  and  equilibrated  with  50  mM  HEPES  (pH  8.0)  and 
500  mM  NaCl,  and  the  column  was  washed  with  2  column 
volumes  of  buffer  after  loading.  The  protein  was  eluted  with  the 
same  buffer  containing  300  mM  imidazole.  Buffer  was 
exchanged  (five  times)  to  50  mM  HEPES  (pH  8.0)  and  150 
mM  NaCl  using  an  ultrafiltration  unit.  Pro¬ 
tein  was  flash-frozen  in  liquid  nitrogen  and  stored  at  —80  °C 
until  it  was  assayed. 

HuGGT  Assays.  Steady  state  kinetic  parameters  were 
measured  in  50  mM  HEPES  (pH  8.0)  and  150  mM  NaCl 
using  L-glutamic  acid  y-(7-amido-4-methylcoumarin)  (l-G1u- 
AMC)  from  Sigma  as  the  substrate.  Rates  were  measured  at 
37  °C  by  monitoring  the  release  of  AMC  continuously 
(excitation  at  380  nm  and  emission  at  460  nm).  The 
fluorescence  intensity  of  AMC  was  proportional  to  the 
concentration  up  to  2  fiM,  and  inner  filter  effect  corrections 
were  not  necessary. 

Cloning  and  Expression  of  K562N,  K562Q,  and  Y403A 
HuGGT  Variants.  WT  HuGGT  and  the  K562N,  K562Q,  and 
Y403A  variants  were  expressed  by  the  same  method  described 
above  using  six-well  plates  and  500  fig  of  plasmid  and  2.5  fiL  of 
lipofectamine  per  well.  Two  controls  were  included:  (l)  cells 
treated  with  lipofectamine  and  no  DNA  and  (2)  untreated  cells. 
The  cells  were  lysed  in  250  fiL  of  lysis  buffer,  and  the  cell 
lysates  were  clarified  by  centrifugation  at  10,000g  for  10  min. 
Protein  concentrations  were  measured  using  a  Bradford  assay 
(Bio-Rad,  Inc.).  Activity  was  measured  using  the  l-G1u-AMC 
substrate  with  or  without  an  acceptor  substrate  and  1:20 
dilutions  of  the  enzyme  (enzyme  diluted  in  assay  buffer). 
Reaction  mixtures  contained  5  fiL  of  the  diluted  enzyme  and 
95  fiL  of  the  substrates  in  assay  buffer  [50  mM  HEPES  (pH  8.0) 
and  150  mM  NaCl].  IC50  values  were  measured  using  the  same 
lysates. 

Modeled  Substrates  and  Inhibitors  of  HuGGT  and 
CapD.  It  is  well-known  that  GGTs  use  glutathione  (GSH),  a 
y-L-Glu-L-Cys-Gly  tripeptide,  as  a  donor  substrate.  However, 
CapD  has  not  been  shown  to  hydrolyze  y-glutamyl-p- 
nitro anilide 9  or  glutathione.3  In  B.  anthracis,  CapD  binds  and 
hydrolyzes  a  polymer  of  y-D-glutamic  acid,  PDGA,  which  forms 
the  capsule.  To  model  PDGA  bound  to  CapD,  we  constructed 
four  derivatives  of  y-D-Glu-D-Glu  dipeptides  in  all  four  capping 
combinations  (i.e.,  acetyl  and/or  N-methyl  capping  groups  at 
the  N-  and  C-termini)  and  refined  their  complexes  using 
MD  simulations  (see  the  Supporting  Information).  Figure  1 
shows  the  glutathione  (GSH)  donor  substrate,  the  HuGGT 
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Figure  1.  Small  molecule  probes  for  GGT  used  in  this  study.  The 
glutathione  HuGGT  substrate  and  B.  anthracis  CapD  substrate  mimic, 
acetylated  di-y-glutamic  acid,  are  shown.  The  OU749  acceptor  site 
inhibitor  of  HuGGT  was  described  by  King  et  al.34 

acceptor  site  inhibitor,  OU749,38  and  the  capped  PDGA  CapD 
substrate  mimic. 

Homology  Modeling  of  HuGGT.  The  atomic  structure  of 
HuGGT  was  built  from  the  crystal  structure  of  EcGGT.  Both 
the  native  structure  and  the  y-glutamyl-bound  structure  of 
EcGGT  were  used  as  templates  (PDB  entries  2DBU  and 
2DBW,  respectively).20  The  restraint-based  comparative 
modeling  program  MODELER  was  used  to  construct  the 
structural  models  of  HuGGT  in  the  native  form  and  y-glutamyl- 
bound  complex.32  Residues  1—34  and  366— 380  were  excluded 
from  the  HuGGT  model,  as  these  are  typically  disordered  in 
the  structures  of  homologues.  A  set  of  20  models  from  random 
generation  of  the  starting  structure  was  minimized,  and  the 
representative  model  was  selected  with  the  lowest  value  of 
MODELER  objective  function.  The  three-dimensional  (3D) 
fold  of  the  generated  models  was  verified  with  PROSA II,40  and 
problematic  regions  were  inspected  via  analysis  of  the  PROSA 
II  energy  profile.  The  resulting  models  were  further  optimized 
with  an  iterative  approach  until  no  significant  improvement  was 
obtained.  The  stereo  chemical  quality  and  protein  structure  of 
the  final  models  were  validated  by  PROCHECK.41 

Molecular  Dynamics  Simulations.  Molecular  dynamics 
(MD)  simulations  were  conducted  for  modeled  systems  in 
explicit  solvent  using  the  AMBER  10  package  and  the 
Parm99SB  force  field.42  The  solvated  protein  systems  were 
subjected  to  energy  minimization  prior  to  MD  simulations  by 
first  minimizing  the  water  molecules  while  holding  the  solute 
frozen  (1000  steps  using  the  steepest  descent  algorithm), 
followed  by  5000  steps  of  conjugate  gradient  minimization  of 
the  whole  system  to  remove  close  contacts  and  to  relax  the 
system.  Bond  lengths  involving  hydrogen  were  constrained  with 
SHAKE,  and  the  time  step  for  all  MD  simulations  was  set  to 
2  fs.  A  nonbonded  cutoff  of  10  A  was  used,  and  the  nonbonded 
pair  list  was  updated  every  25  steps.  Periodic  boundary 
conditions  were  applied  in  simulations  using  the  particle  mesh 
Ewald  (PME)  method  to  calculate  the  long-range  electrostatic 
interactions.  The  simulated  system  was  first  subjected  to  a 
gradual  temperature  increase  from  0  to  300  K  over  100  ps  and 
then  equilibrated  for  500  ps  at  300  K,  followed  by  production 
runs  with  a  total  length  of  2  ns.  Constant  temperature  and 
pressure  (300  K  and  1  atm,  respectively)  were  maintained  using 
the  Berendsen  coupling  algorithm  with  a  time  constant  for 
heat-bath  coupling  of  0.2  ps.  The  resulting  trajectories  were 
analyzed  using  the  PTRAJ  module  from  AMBER  The  root- 
mean-square  deviations  (rmsds)  of  the  backbone  were 
calculated  from  the  trajectories  at  1  ps  intervals,  with  the 


initial  structure  as  the  reference.  The  root-mean-square 
positional  fluctuations  (rmsfs)  were  calculated  on  the  basis  of 
superposition  of  all  C„  atoms  using  the  averaged  structure  of 
trajectories  as  the  reference  and  compared  to  the  experimental 
B  factors. 

Ensemble  Docking  and  Binding  Free  Energy  Calcu¬ 
lations.  The  AutoDock-based  DOVIS  program43  was  used  for 
docking  the  substrates  to  the  modeled  structure  of  HuGGT  and 
CapD.  Multiple  conformations  associated  with  the  active  site 
residues  were  obtained  from  MD  simulations  and  used  for  all 
docking  studies.  The  active  site  of  the  protein  was  defined  by  a 
grid  of  70  X  70  X  70  points  with  a  grid  spacing  of  0.375  A 
centered  at  the  catalytic  residue  Thr.  The  Lamarckian  genetic 
algorithm  (LGA)  in  AutoDock  was  applied  with  50  runs,  and 
the  maximal  number  of  energy  evaluations  was  set  to  2  X  106. 
Results  differing  by  <  1.0  A  in  the  positional  rmsd  of  the 
substrate  were  clustered,  and  the  one  with  the  most  favorable 
free  energy  of  binding  represented  the  final  binding 
conformations.  The  optimal  binding  complexes  were  subjected 
to  a  stepwise  energy  minimization  and  MD  simulations  in 
explicit  solvent  as  described  above.  The  binding  free  energies 
were  calculated  using  the  MM-GBSA  method  as  described 
previously.44  A  set  of  300  snapshots  was  extracted  from 
trajectories  of  binding  complexes  at  10  ps  intervals  from  the  last 
3  ns  of  each  MD  simulation.  The  molecular  mechanics  inter¬ 
action  energies  were  calculated  with  the  SANDER  module  in 
the  Amber  10  package.  The  entropy  was  estimated  through 
normal-mode  analysis  using  the  NMODE  module. 

■  RESULTS 

Kinetic  Analysis  of  the  Acceptor  Substrate  Prefer¬ 
ences  in  CapD.  Using  the  purified  poly-D-glutamate  (PDGA) 
capsule  from  the  Ames  strain  of  B.  anthracis,  we  initially 
observed  a  faster  rate  of  B.  anthracis  PDGA  capsule  degradation 
in  the  presence  of  cell  culture  medium  versus  buffer.  Analysis  of 
PDGA  capsule  degradation  by  agarose  gel  electrophoresis 
indicated  that  the  transpeptidase  activity  was  greater  than  the 
hydrolase  activity.  The  components  of  the  media  were  tested 
individually;  some  but  not  all  amino  acids  were  able  to  enhance 
capsule  cleavage  (Figure  SI  of  the  Supporting  Information).  To 
quantify  the  rate  of  transpeptidation  and  to  identify  any  unique 
amino  acid  acceptor  substrate  preferences  of  the  bacterial  GGT, 
CapD,  that  differ  from  those  of  rat  GGT  (and  presumably 
HuGGT),  the  steady  state  kinetic  parameters  were  measured 
for  each  amino  acid  by  globally  fitting  three-dimensional  data 
sets  (Table  l).  Both  FRET  donor  and  acceptor  substrate 
concentrations  were  varied  and  yielded  a  set  of  parallel  lines  in 
a  double-reciprocal  plot  consistent  with  a  ping-pong  bi-bi 
mechanism  (binary  complex)  (Figure  S2  of  the  Supporting 
Information). 

From  the  three-dimensional  kinetic  analysis,  we  found  that 
the  Kms  app  values  for  the  FRET  donor  substrate  varied  for  each 
of  the  amino  acid  acceptor  substrates.  For  a  ping-pong 
mechanism,  the  binding  of  the  donor  substrate  and  formation 
of  the  covalent  intermediate  are  expected  to  be  the  same 
regardless  of  the  acceptor  substrate.  Thus,  the  variation  likely 
arises  from  competitive  inhibition  by  the  acceptor  substrate; 
this  was  previously  suggested  by  Thompson  and  Meister  in 
their  studies  of  rat  GGT.26 

Comparison  of  the  bimolecular  rate  constants,  fccal/JCma,  for 
each  amino  acid  acceptor  substrate  showed  the  order  of  the 
CapD  acceptor  substrate  preference:  L-Cys  >  l-G1u  >  L-Gln  > 
d-G1u  =  L-Met  (Table  l).  Of  the  acceptors  tested,  L-Cys  yielded 
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Table  1.  Acceptor  Substrate  Preferences  for  CapD-MBP  in  25  mM  HEPES  (pEI  7.4 )  and  0.1%  Tween  20  at  22  +  3  °C“ 


acceptor  substrate 

jFCmsapp  (FRET  donor  substrate)  (//M) 

Rma  (acceptor)  (//M) 

Vmax  (Units/mg) 

fccat  (min-1) 

fccat/Rma  (mm  1  //M-1) 

h2o 

4  ±  2 

- 

0.0019  ±  0.0003 

0.19  ±  0.03 

- 

L-alanine 

6  +  2 

5000  ±  1000 

0.011  ±  0.002 

1.1  ±  0.2 

0.00022  ±  0.00006 

L-arginine 

50  ±  10 

40000  ±  10000 

0.07  ±  0.02 

7  ±  2 

0.00018  ±  0.00006 

L-asparagine 

9  ±  3 

15000  ±  5000 

0.035  ±  0.008 

3.5  ±  0.8 

0.00023  ±  0.00009 

L-aspartate 

7  ±  3 

300  ±  100 

0.0044  ±  0.0008 

0.44  ±  0.08 

0.0015  ±  0.0006 

L-cysteine 

12  ±  2 

170  ±  20 

0.023  ±  0.002 

2.3  ±  0.2 

0.014  ±  0.002 

D-glutamate 

8  ±  3 

150  ±  40 

0.006  ±  0.001 

0.6  ±  0.1 

0.004  ±  0.001 

L-glutamate 

22  ±  7 

400  ±  100 

0.028  ±  0.006 

2.8  ±  0.6 

0.007  ±  0.002 

L-glutamine 

14  ±  9 

400  ±  300 

0.02  ±  0.01 

2  ±  1 

0.005  ±  0.005 

glycine 

13  ±  11 

5000  ±  4000 

0.008  ±  0.004 

0.8  ±  0.4 

0.0002  ±  0.0002 

L-histidine 

3  ±  2 

210  ±  150 

0.006  ±  0.001 

0.6  ±  0.1 

0.003  ±  0.002 

L-isoleucine 

7  ±  2 

25000  ±  6000 

0.019  ±  0.003 

1.9  ±  0.3 

0.00008  ±  0.00002 

L-leucine 

8  ±  4 

600  ±  200 

0.016  ±  0.004 

1.6  ±  0.4 

0.003  ±  0.001 

L-lysine 

19  ±  9 

800  ±  500 

0.010  ±  0.003 

1.0  ±  0.3 

0.0013  ±  0.0009 

L-methionine 

2  ±  1 

80  ±  30 

0.0029  ±  0.0005 

0.3  ±  0.05 

0.004  ±  0.002 

L-phenylalanine 

8  ±  4 

1500  ±  500 

0.015  ±  0.004 

1.5  ±  0.4 

0.0010  ±  0.0004 

L-proline 

NAb 

L-serine 

7.8  ±  0.7 

990  ±  80 

0.021  ±  0.001 

2.1  ±  0.1 

0.0021  ±  0.0002 

L-threonine 

6  ±  4 

1400  ±  900 

0.010  ±  0.003 

1.0  ±  0.3 

0.0007  ±  0.0005 

L-tryptophan 

13  ±  6 

1500  ±  500 

0.013  ±  0.003 

1.3  ±  0.2 

0.0009  ±  0.0003 

L-valine 

20  ±  10 

32000  ±  17000 

0.03  ±  0.01 

3.0  ±  1.0 

0.00009  ±  0.00006 

aL-Tyrosine  was  not  tested  because  of  its  poor  solubility  at  neutral  pH.  Steady  state  kinetic  parameters  for  CapD  without  MBP  were  measured  using 
L-Cys  [fCms  app  (substrate)  =  5  ±  2  /^M,  Kma  =  220  ±  60  fdM,  Vmax  =  0.024  ±  0.004  Units/mg,  and  fccat  =  1.3  ±  0.2  min  *]  and  were  found  to  be 
comparable  to  those  measured  for  CapD-MBP.  bNo  significant  increase  in  activity  was  observed  with  50  mM  L-proline.  L-Proline  does  not  have  a 
primary  amine  that  can  function  as  an  acceptor  in  the  reaction. 


Amino  Acid  Acceptor 

Figure  2.  Comparison  of  acceptor  substrate  preferences  for  rat  GGT 
(published  values  from  ref  26)  and  B.  anthracis  CapD  (from  Table  l). 
Both  GGT  enzymes  prefer  L-Cys  as  their  acceptor;  however,  CapD, 
unlike  rat  GGT,  can  also  utilize  l-G1u  and  other  amino  acids  as 
acceptor  substrates. 


the  highest  fccat/JCma  and  L-Met  had  the  lowest  JCma.  CapD, 
therefore,  had  an  acceptor  preference  similar  to  that  of  rat  GGT 
(which  is  79%  identical  to  the  HuGGT  sequence)  and  has  been 
reported  to  have  the  same  acceptor  substrate  preference  for 
L-Cys  and  the  highest  affinity  for  L-Met  (Figure  2). 13,26  Rat 
GGT  prefers  the  L-Cys  acceptor  substrate  most  likely  because  of 
its  similarity  to  its  donor  substrate  (glutathione)  and  product 
(cysteinyl-glycine).  However,  for  CapD,  this  result  was  unexpected 
because  glutathione  has  not  been  shown  to  be  a  donor  substrate.9 
The  common  acceptor  preference  observed  between  CapD  and 
rat  GGT  suggests  that,  although  their  donor  substrate  preferences 
are  distinctly  different,  CapD  and  rat  GGT  may  share  similar 
binding  interactions  within  their  acceptor  sites. 


Further  examination  of  the  CapD  acceptor  substrate 
preferences  showed  that  CapD  could  utilize  a  broader  range 
of  acceptor  substrates  than  rat  GGT  (Figure  2).  Because  CapD 
hydrolyzes  an  unusual  poly-D-glutamate  donor  substrate,3  we 
tested  the  negatively  charged  d-G1u  amino  acid  as  an  acceptor 
substrate  in  transpeptidation.  A  2.7-fold  lower  Rma  was  found 
for  d-G1u  compared  to  that  for  l-G1u.  However,  the  turnover 
number,  fcca0  with  d-G1u  was  4.7-fold  lower  than  that  with 
l-G1u.  Overall,  the  specificity  constant,  fccat/Rma)  for  l-G1u  was 
1.8-fold  higher  than  for  d-G1u  (Table  l),  indicating  that  the 
acceptor  site  of  CapD  does  not  strongly  prefer  D-amino  acids 
over  L-amino  acids.  Rat  GGT  and  EcGGT  have  strict 
stereospecificity  for  L-amino  acid  acceptors  over  D-amino 
acids.  The  lack  of  stereospecificity  in  CapD  demonstrates  that 
the  acceptor  site  of  CapD,  while  similar  in  some  respects,  is 
distinctly  different  from  that  of  HuGGT.  To  understand  the 
lack  of  stereospecificity  and  the  acceptor  substrate  preference 
for  negatively  charged  amino  acids,  we  utilized  both  site- 
directed  mutagenesis  and  molecular  modeling  of  the  d-  and 
L-amino  acid  acceptor  substrates. 

Site-Directed  Mutagenesis  Studies  of  the  CapD 
Acceptor  Site.  Examination  of  the  putative  acceptor  substrate 
site  in  CapD  suggested  that  Arg-432  and  Arg-520  may  both  be 
involved  in  binding  the  acceptor  substrate  (Figure  3A  and 
Figures  S3  and  S4  of  the  Supporting  Information).  For  CapD, 
the  L-Arg  acceptor  substrate  yielded  the  highest  fccat.  The  higher 
turnover  number  may  be  due  to  faster  product  release,  as  the 
ICma  of  this  acceptor  [40  ±  10  mM  (Table  l)]  was  the  highest 
of  the  amino  acids  tested.  Docking  indicated  that  the  poor 
binding  of  this  acceptor  may  be  due  to  Arg-432,  which  is 
predicted  to  be  within  the  putative  acceptor  site  (Figure  3A). 
This  result  was  similarly  observed  for  EcGGT,19  suggesting  that 
the  two  bacterial  GGTs  share  a  similar  residue  within  their 
acceptor  substrate  sites  (Figure  3B).  EcGGT  also  has  an 
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Figure  3.  (A)  Predicted  binding  models  of  preferred  amino  acid  acceptors  of  CapD.  (B)  Structural  and  sequence  alignment  of  the  putative  acceptor 
site  of  CapD  with  those  of  other  GGTs.  The  numbering  at  the  top  corresponds  to  that  of  CapD  (purple);  the  numbering  at  the  bottom  of  the 
alignment  corresponds  to  that  of  the  human  enzyme  (yellow).  The  residues  corresponding  to  EcGGT  are  colored  blue. 


arginine  at  the  position  homologous  to  Arg-432  in  CapD,  while 
in  HuGGT,  a  neutral  glutamine  is  present.  We  hypothesized 
that  differences  between  the  bacterial  and  human  acceptor 
substrate  preferences  could  arise  from  this  charged  residue. 

To  further  validate  our  predicted  acceptor  site,  we  sought 
experimental  evidence  and  characterized  three  CapD  variants, 
N373A,  R432A,  and  R520S  (Table  2).  The  Asn-373  variant  was 
selected  because  it  is  approximately  180°  from  Arg-432. 
Analysis  of  the  kinetic  parameters  of  each  of  the  variants 
allowed  the  orientation  of  the  side  chain  of  the  acceptor  amino 
acid  in  our  models  to  be  verified.  To  establish  the  orientation  of 
the  COO-  group  of  the  acceptor  substrate,  which  could  point 
toward  either  Arg-520  or  Arg-432,  we  also  selected  the  R520S 
variant.  In  HuGGT,  the  residue  homologous  to  Arg-520  is 
conservatively  substituted  with  Lys-562.  Previously,  no  role  for 
this  positively  charged  residue  had  been  established  in  catalysis. 

Site-directed  mutagenesis  of  Arg-520  to  serine  showed 
a  >  95-fold  decrease  in  fccat  (hydrolysis)  consistent  with  a  role 


for  Arg-520  in  catalysis  and/or  product  release  (Table  2).  The 
activity  of  the  R520S  variant  was  not  stimulated  by  the  addition 
of  an  amino  acid  acceptor  substrate,  and  no  detectable  activity 
for  this  variant  could  be  measured  (a  1400-fold  decrease  in  fccat 
compared  to  that  of  WT  CapD  using  L-Cys  as  the  acceptor 
substrate).  Thus,  the  role  of  Arg-520  may  be  to  interact  with 
the  carboxylate  of  the  amino  acid  acceptor  substrate  and  func¬ 
tion  as  an  electrostatic  catalyst. 

Examination  of  a  second  variant,  R432A,  showed  a  smaller 
6-fold  decrease  in  fccat  and  no  significant  change  in  the  apparent 
Kms  of  the  donor  substrate  when  water  functioned  as  the 
acceptor  substrate.  The  4.9-fold  decrease  in  fccat/ Km  for  the 
hydrolase  activity  suggests  a  minor  role  for  Arg-432  in  catalysis. 
However,  like  the  R520S  variant,  no  significant  stimulation  of 
the  activity  was  observed  when  the  Gly,  L-Asp,  d-G1u,  or  l-G1u 
acceptor  was  present  at  5  mM.  Rather,  a  decrease  in  activity  was 
observed  when  d-G1u,  l-G1u,  or  L-Asp  was  added.  Further 
examination  of  the  R432A  variant  showed  that  l-G1u  could  now 
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Table  2.  Steady  State  Kinetic  Parameters  for  CapD  Acceptor  Site  Mutants 


enzyme 

acceptor  substrate 

JCmsapp  (FRET  donor  substrate)  (//M) 

Kma  (acceptor)  (//M) 

Vmax  (Units/mg) 

fccat  (min-1) 

fccat/.fCma  (min  1  //IVT1 

R432A 

h2o 

3.2  ±  0.7 

0.00031  ±  0.00002 

0.031  ±  0.002 

h2o 

—a 

0.00031  ±  0.00001 

0.031  ±  0.001 

glycine 

10  ±  10 

>20000 

>0.0010 

>0.10 

L-arginine 

10  ±  6 

1500  ±  1400 

0.0024  ±  0.0008 

0.24  ±  0.08 

0.0002  ±  0.0002 

L-aspartate 

—a 

>5000 

>0.00007 

>0.007 

L-cysteine 

—a 

>5000 

>0.0005 

>0.05 

D-glutamate 

—a 

>5000 

>0.00007 

>0.007 

L-glutamate 

—a 

>5000 

>0.0001 

>0.01 

N373A 

h2o 

3  ±  1 

0.00050  ±  0.00008 

0.050  ±  0.008 

L-cysteine 

10  ±  3 

500  ±  100 

0.008  ±  0.001 

0.8  ±  0.1 

0.0016  ±  0.0004 

D-glutamate 

5  ±  2 

60  ±  40 

0.006  ±  0.002 

0.6  ±  0.2 

0.010  ±  0.007 

L-glutamate 

15  ±  8 

120  ±  70 

0.015  ±  0.006 

1.5  ±  0.6 

0.013  ±  0.009 

R520S 

h2o 

_b 

<0.00002 

L-glutamate 

—a 

>5000 

<0.00002 

“Activity  was  measured  using  40  j.iM  FRET  substrate  and  <5  mM  acceptor  amino  acid.  bNo  significant  activity  detected.  Limit  of  detection  reported. 


act  as  a  competitive  inhibitor  with  a  K,  of  0.5  ±  0.3  mM.  The 
JCma  for  L-Arg  with  the  R432A  variant  decreased  dramatically 
from  40  to  1.5  mM  [26-fold  (Table  2)],  corroborating  the 
orientation  of  the  acceptor  substrate  side  chain  toward  R432. 
The  fccat  correspondingly  decreased  29-fold,  consistent  with  a 
slower  rate  of  product  release.  Thus,  the  role  of  Arg-432 
appears  to  be  primarily  in  binding  the  R  group  of  the  L-amino 
acid  acceptor  substrates. 

Conversely,  D-amino  acid  acceptor  substrates  appear  to 
employ  Arg-432  in  binding  the  a-carboxylate,  as  opposed  to 
Arg-520.  The  two  alternate  binding  modes  for  d-  and  l- 
acceptor  substrates  are  shown  in  Figure  3.  The  two  arginines 
could  account  for  the  lack  of  stereospecificity  in  CapD.  EcGGT 
and  rat  GGT  have  strict  stereospecificity  for  L-acceptor 
substrates  and  do  not  have  both  arginines  (Figure  3B). 
Consistent  with  our  models,  the  R432A  variant  was  unable  to 
utilize  d-G1u  or  l-G1u  as  an  acceptor  substrate  (Table  2).  The 
sulfur  atom  of  Met  is  known  to  form  hydrogen  bonds  with  NFI 
groups.45-47  A  hydrogen  bond  between  the  Met  sulfur  and  Arg- 
432  side  chain  may  in  part  account  for  the  tight  binding  of  this 
acceptor  (Figure  3). 

A  third  variant,  N373A,  was  examined.  Asn-373  was  also 
predicted  to  be  near  the  putative  acceptor  substrate  site  (Figure  3). 
However,  the  N373A  variant  showed  an  only  2.9-fold  decrease  in 
fccat  when  water  functioned  as  the  acceptor.  The  JCma  of  L-Cys  was 
also  2.9-fold  higher  than  the  Kmi  of  the  L-Cys  measured  for  WT 
CapD  and  was  not  significantly  different,  suggesting  that  the 
acceptor  side  chain  is  not  oriented  toward  this  residue.  The 
measured  Kma  values  for  d-G1u  and  l-G1u  were  also  comparable  to 
those  of  the  wild  type,  suggesting  that  the  Asn-373  side  chain  does 
not  play  a  significant  role  in  binding  either  stereoisomer. 

Structural  Properties  of  Human  GGT  and  CapD.  To 
improve  our  understanding  of  the  structural  basis  of  the 
acceptor  preference  associated  with  HuGGT  and  CapD,  we 
applied  computational  approaches  to  probe  the  binding 
specificities  of  the  substrate  and  inhibitors  in  the  donor  and 
acceptor  sites.  Because  the  structure  of  HuGGT  is  not  available, 
we  modeled  it  on  the  basis  of  homologous  bacterial  GGT.  In 
addition  to  HuGGT,  the  structure  of  rat  GGT  was  also 
modeled,  which  possesses  an  active  site  and  ternary  structure 
identical  to  those  of  HuGGT  (Figure  S5  of  the  Supporting 
Information).  HuGGT  shares  the  highest  level  of  sequence 
identity  (~31%)  with  EcGGT  and  levels  of  20—28%  with  other 
bacterial  GGTs  from  H.  pylori,  Bacillus  subtilis  (BsGGT),  and 


B.  anthracis.  The  modeled  structure  of  HuGGT  showed  an  overall 
rmsd  for  Ca  atoms  of  0.36  A  with  respect  to  the  homologue, 
EcGGT,  and  values  of  0.83,  1.19,  and  1.29  A  compared  to 
the  experimental  structures  of  HpGGT,  BsGGT,  and  CapD, 
respectively. 

Structural  analysis  revealed  that  the  same  3D  fold  is  shared 
between  human  and  bacterial  GGTs  (Figure  4).  The  most 
varied  regions  are  located  on  the  large  40  kDa  subunit  and 
surface  loops  at  the  terminal  end,  while  the  catalytic  binding 
site  is  mainly  formed  by  the  small  20  kDa  subunit,  including  the 
catalytic  residue  Thr-381  and  the  conserved  oxyanion  hole 
motif  of  Gly-483  and  Gly-484  in  HuGGT.  Compared  to  CapD, 
the  most  distinct  structural  difference  is  the  fact  that  HuGGT 
contains  a  substrate  binding  loop  at  the  donor  site  (residues 
427—438),  whereas  such  a  loop  is  absent  in  B.  anthracis  CapD 
as  well  the  GGT  in  B.  subtilis  and  S.  aureus 48  Interestingly, 
another  B.  subtilis  GGT  containing  a  donor  substrate  binding 
loop  has  also  been  identified.49  The  substrate-binding  loop  is 
found  in  many  other  bacterial  GGTs  and  plays  an  important 
role  in  substrate  recognition.2-  The  lack  of  the  loop  in  CapD  is 
strongly  associated  with  the  substrate  and  inhibitor  binding 
specificities,  as  shown  in  the  following  modeling  studies. 
Another  notable  difference  between  HuGGT  and  bacterial 
GGTs  is  that  HuGGT  possesses  a  “gating  residue”  Phe-433 
within  the  substrate  binding  loop,  while  it  is  typically  a  tyrosine 
in  bacterial  GGTs. 

Closer  examination  of  the  HuGGT  and  CapD  structures 
suggested  that  the  putative  acceptor  binding  site  is  located  in 
the  region  along  the  groove  of  the  donor  binding  site  (Figure  4). 
A  large  pocket  is  found  in  the  proximity  of  the  catalytic  resi¬ 
due  and  y-glutamyl  intermediate  that  possibly  accommo¬ 
dates  various  acceptor  amino  acids  for  the  transpeptidation 
reaction.  In  contrast  to  the  conserved  donor  site,  the  acceptor 
site  residues  varied  greatly  between  human  and  bacterial  GGTs. 
The  binding  pocket  in  HuGGT  is  formed  by  a  number  of 
polar  residues  such  as  Asp-46,  Asn-79,  His-81,  Ser-82,  Gln-476, 
Tyr-403,  and  Lys-562.  The  residues  equivalent  to  Gln-476  and 
Lys-562  are  substituted  with  Arg-432  and  Arg-520,  respectively, 
in  CapD;  these  residues  were  shown  to  be  involved  in  acceptor 
binding  in  our  site-directed  mutagenesis  study.  Residues  Arg- 
432  and  Gln-476  are  located  on  top  of  the  donor-binding 
pocket  and  are  predicted  to  form  hydrogen  bond  interactions 
with  the  acceptor  substrate  (Figure  3).  On  the  other  hand, 
Arg-520  and  Lys-562  point  toward  the  acceptor  site  from  the 
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Figure  4.  (A)  Structural  model  of  HuGGT  (magenta  ribbon  and  red  loop)  in  comparison  with  CapD  (blue  ribbon,  PDB  entry  3G9K).25  (B) 
HuGGT  illustrated  in  surface  representation  (white).  The  HuGGT  substrate-binding  loop  (red),  nucleophilic  residue  T381  (magenta),  and  docked 
glutamyl  group  (green)  are  shown  as  sticks.  The  CapD  structure  is  shown  as  ribbons  (blue). 


C-terminal  loop  of  the  small  subunit.  This  residue  appears  to  be 
highly  flexible  as  two  conformations  were  found  in  the  X-ray 
crystal  structure  of  CapD,25  and  various  conformations  were 
seen  in  other  bacterial  GGTs.  The  motion  of  the  C-terminal 
loop  in  HpGGT  and  the  functional  role  in  substrate  bind¬ 
ing  have  been  reported  recently.23  Notably,  both  HuGGT  and 
CapD  possess  an  aromatic  residue  (Phe-374  and  Tyr-403, 
respectively)  located  in  the  front  of  the  acceptor  binding 
pocket,  while  a  smaller  residue  is  found  in  EcGGT  (Thr-413) 
and  HpGGT  (Ser-402)  with  a  rather  open  putative  acceptor 
site.  It  is  likely  that  this  residue  plays  a  functional  role  similar  to 
that  of  the  gating  residue  in  the  donor  substrate-binding  loop, 
shielding  the  acceptor  substrate  or  inhibitor  from  solvent. 

Dynamics  of  the  Acceptor  Substrate  Site  Occur  in 
Cooperation  with  the  Donor  Substrate  Site.  To 
investigate  the  dynamic  properties  of  HuGGT  and  CapD,  we 
performed  MD  simulations  for  the  apo  and  y-glutamyl-enzyme 
intermediate  (acylated)  complexes  in  explicit  solvent.  Figure  5 
depicts  the  atomic  positional  fluctuations  of  residues  (Ca) 
for  both  systems  calculated  from  2  ns  of  MD  simulations. 
The  calculated  B  factors  were  generally  in  agreement  with  the 
experimentally  measured  data  from  the  glutamyl-enzyme 
intermediate  of  EcGGT.20  Larger  atomic  fluctuations  were 
mainly  observed  in  loop  regions  on  the  protein  surface  and 
termini  of  two  subunits  (except  for  the  N-terminal  Thr-381  in 
HuGGT  and  Thr-352  in  CapD).  Both  the  donor  and  acceptor 
sites  of  HuGGT  and  CapD  exhibited  a  higher  level  of 
conformational  flexibility  in  the  apo  form,  and  the  magnitude 
of  the  fluctuations  decreased  in  the  glutamyl-enzyme 
intermediate  complex.  The  putative  acceptor  site  appears  to 
be  more  flexible  than  the  donor  site.  Analysis  showed  that 
a  well-formed  hydrogen  bonding  network  rigidifies  the 
donor  substrate  binding  site,  while  the  acceptor  site  residues 
remain  flexible.  Interestingly,  the  presence  of  the  y-glutamyl 
moiety  did  not  affect  the  CapD  acceptor  substrate  binding  site 
significantly,  implying  that  the  effect  of  the  donor  substrate 
binding  on  the  acceptor  site  is  different  from  that  of  HuGGT. 
It  is  worth  noting  that  two  small  loops  of  HuGGT  at  the 
acceptor  site,  L2  (residues  459—464)  and  L3  (residues  549— 
553),  were  found  to  exhibit  large  motions.  In  CapD,  helix  a5 
(residues  167—179)  in  the  back  of  the  acceptor  binding  site  is 
attenuated  upon  acylation  of  the  enzyme. 

Substrate  binding  loop  LI  (residues  423—438)  of  HuGGT 
exhibited  large  fluctuations  in  the  apo  state  as  expected.  The 
result  is  consistent  with  that  observed  in  the  apo  structure  of 


Figure  5.  B  factors  of  C„  atoms  computed  from  MD  simulations  of 
HuGGT  (A)  and  CapD  (B).  Results  for  both  apo  (blue  solid  line)  and 
acylated  (red  dashed  line)  forms  are  presented.  Yellow  and  green 
horizontal  bars  atop  each  plot  show  the  boundaries  of  the  large  (L) 
and  small  (S)  domains.  In  the  acylated  enzymes,  the  B  factors  of  the 
acyl  C„  atoms  precede  those  of  the  catalytic  threonine  residues  (T381 
in  HuGGT  and  T352  in  CapD).  The  regions  of  high  flexibility  in 
CapD  that  are  labeled  with  red  bars  correspond  to  modeled  loops  the 
density  of  which  is  not  defined  in  the  X-ray  crystallographic  structure 
(PDB  entry  3G9K). 


EcGGT  and  HpGGT.  However,  the  y-glutamyl  moiety  did  not 
significantly  reduce  the  flexibility  of  the  LI  loop  in  HuGGT. 
Analysis  showed  that  the  gating  residue,  Tyr-444  in  EcGGT 
and  Tyr-433  in  HpGGT,  within  the  loop  forms  a  hydrogen 
bond  with  the  amide  of  the  glutamyl  moiety,  whereas  HuGGT 
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Figure  6.  (A)  Binding  interactions  of  GSH  in  the  active  site  of  HuGGT.  Residues  involved  in  protein— ligand  interactions  are  shown  as  sticks  (cyan 
and  red).  For  the  sake  of  clarity,  hydrogen  atoms  of  the  ligand  have  been  omitted.  Pairs  of  atoms  involved  in  hydrogen  bonding  are  connected  with 
dashed  lines.  (B)  Relative  activity  of  WT  HuGGT  and  the  K562N,  K562Q,  and  Y403A  variants  with  or  without  40  mM  glycyl-glycine.  (C)  No 
increase  in  activity  was  observed  with  increasing  concentrations  of  the  L-Cys  acceptor  substrate  for  the  K562N,  K562Q,  and  Y403A  variants.  (D) 
Binding  interactions  of  acetylated  di-y-glutamic  acid  (green)  with  CapD  (magenta). 


containing  a  Phe-433  gating  residue  lacks  such  a  stabilizing 
binding  interaction.  As  a  result,  the  substrate  binding  loop  of 
HuGGT  in  the  substrate-bound  complex  exhibited  larger 
fluctuations  than  EcGGT  and  HpGGT.  The  conformational 
change  of  the  substrate  binding  loop,  from  open  to  closed,  is 
thought  to  facilitate  substrate  recognition  and  catalysis. 
Therefore,  the  dynamics  and  higher  flexibility  of  the  substrate 
binding  loop  of  HuGGT  are  likely  associated  with  its  higher 
transpeptidase  activity  (faster  product  release). 

Differences  in  Donor  Substrate  Binding  Specificities 
of  HuGGT  and  CapD.  GGTs  that  lack  the  donor  substrate 
binding  loop  have  not  been  shown  to  utilize  glutathione  as  a 
substrate.  The  absence  of  this  loop  suggests  a  loss  of  key 
interactions  with  the  GSH  substrate  of  HuGGT  and  its 
glutamine  analogue  inhibitors  such  as  acivicin,  which  inhibits 
CapD  poorly  (IC50  >  20  mM).  To  understand  why  GSH  does 
not  function  as  a  substrate  of  CapD,  we  investigated  the 
binding  specificity  of  GSH  with  HuGGT.  Figure  6A  shows  the 
predicted  complex  of  GSH  bound  in  the  active  site  of  HuGGT. 
As  expected,  GSH  adopted  a  binding  mode  similar  to  that 
observed  in  the  experimental  structure  of  the  HpGGT  T380A 
variant  in  complex  with  S-(nitrobenzyl)glutathione.22  Impor¬ 
tantly,  GSH  interacts  with  both  donor  and  acceptor  binding 
sites  of  HuGGT.  The  '/-glutamyl  moiety  of  GSH  was  buried 
deep  in  the  pocket  of  the  donor  substrate  site  and  formed  an 
extensive  hydrogen  bonding  network  with  residues  Arg-103, 
Asn-401,  Glu-419,  Asp-423,  Ser-451,  and  Ser-452.  These 
residues  are  highly  conserved  among  GGTs  that  contain  the 
donor  substrate  binding  loop,  implying  their  primary  functional 
role  in  donor  substrate  recognition.  Moreover,  reported 


mutational  studies  of  HuGGT  have  supported  the  involvement 
of  Asp-423,  Ser-451,  and  Ser-452  in  donor  substrate  bind¬ 
ing.29,50  Asp-423  in  HuGGT  formed  H-bond  and  electrostatic 
interactions  with  the  a-amino  group  of  the  /-glutamyl  moiety. 
In  contrast,  Asp-423  in  CapD  is  substituted  with  an  asparagine. 
The  loss  of  this  important  interaction,  together  with  the  lack  of 
a  substrate  binding  loop,  likely  explains  why  GSH  does  not 
stably  bind  within  the  donor  substrate  site  or  function  as  a 
substrate  in  CapD. 

Further  analysis  of  the  GSH  binding  mode  indicated  that  the 
Cys-Gly  leaving  group  pointed  toward  the  putative  acceptor  site 
of  HuGGT  and  formed  two  interesting  interactions  (Figure  6A). 
One  potential  hydrogen  bond  was  found  between  the  side 
chain  of  the  Cys  acceptor  substrate  and  Asn-476,  similar  to  the 
predicted  binding  model  of  acceptor  L-Cys  bound  in  CapD,  as 
well  as  the  structure  of  the  HpGGT  T380A  variant  bound  to  a 
glutathione  analogue.22  Another  interesting  interaction  was 
observed  with  residues  Lys-562  and  Tyr-403,  which  formed 
hydrogen  bonds  and  strong  charge— charge  electrostatic 
interactions  with  the  glycyl  moiety  of  GSH.  In  contrast, 
residues  Thr-413  and  Ser-572  occupied  these  positions  in 
EcGGT  and  did  not  participate  in  any  interactions  with  the 
glycyl  moiety  in  the  modeled  GSH— EcGGT  complex.  MD 
simulations  showed  that  the  Cys-Gly  leaving  group  of  GSH  in 
HuGGT  was  stabilized  because  of  well-formed  interactions  with 
residues  Lys-562  and  Tyr-403.  This  may  benefit  the 
nucleophilic  attack  in  the  transpeptidation  reaction  involving 
glycyl  dipeptides  and  also  account  for  the  higher  transpeptidase 
activity16  compared  to  that  of  EcGGT.19 
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CapD,  unlike  other  GGTs,  utilizes  a  polymer  of  y-D-glutamic 
acid  as  a  substrate.  To  probe  the  substrate  binding  specificities, 
we  modeled  four  derivatives  of  y-D-Glu-D-Glu  dipeptides  in  all 
four  capping  combinations  (i.e.,  acetyl  and/or  N-methyl 
capping  groups  at  N-  and  C-termini,  respectively)  and  docked 
them  into  the  active  site  of  CapD.  Figure  6  shows  that  PDGA 
could  interact  with  both  the  donor  and  acceptor  substrate  sites 
and  bind  between  Arg-113  in  the  donor  substrate  site  and  Arg- 
432  and  Arg-520  in  the  putative  acceptor  substrate  site.  The 
substrate  can  bind  in  two  different  ways.  Specifically,  the 
substrate  chain  can  swap  its  N-  and  C-termini.  However,  the 
conformer  in  which  the  Arg-113  side  chain  is  closer  to  the 
N-terminus  appears  to  correspond  to  the  productive  confor¬ 
mation,  as  it  satisfies  the  oxyanion  hole  constraint  and  has  a 
favorable  distance  to  the  Or  atom  of  the  nucleophilic  Thr-352 
residue. 

Site-Directed  Mutagenesis  Studies  of  the  HuGGT 
Acceptor  Site.  To  obtain  experimental  evidence  for  the 
predicted  residues  involved  in  substrate  and  inhibitor  binding, 
we  conducted  site-directed  mutagenesis  studies  of  the  HuGGT 
acceptor  site.  The  His-tagged  WT  HuGGT  was  purified  and 
assayed  using  the  l-G1u-AMC  substrate.  For  the  hydrolysis 
reaction,  a  Kms  of  1.7  ±  0.2  f,iM  was  measured  and  agreed  well 
with  the  Kms  previously  measured  by  Han  et  al.18  (fCms  =  4  //M) 
at  pH  8.0  (25  °C)  using  the  same  substrate  in  the  absence  of  an 
acceptor  substrate.  For  the  transpeptidase  reaction,  Kms  =  16  ± 
9  /(M  for  the  l-G1u-AMC  donor  substrate  and  K ma  =  2200  ± 
1500  f.i M  for  the  L-Cys  acceptor  substrate.  The  Km3  of  L-Cys  for 
HuGGT  was  6-fold  higher  than  the  ICm3  of  L-Cys  reported  for 
rat  GGT  (_Kma  =  340  ±  60  ^M)26  at  pH  8.0  and  37  °C.  The 
slightly  higher  Km3  of  L-Cys  may  be  due  to  Ala-472,  which  is  a 
serine  in  rat,  pig,  and  mouse  GGT. 

To  analyze  the  effects  of  the  mutations,  WT  HuGGT  and  the 
K562N,  K562Q,  and  Y403A  variants  were  expressed  in 
quadruplicate  in  six-well  plates.  Lysates  of  the  COS-1  cells 
were  assayed  with  the  l-G1u-AMC  donor  substrate  (Figure  6). 
Two  controls  were  included:  (l)  no  DNA  and  lipofectamine 
and  (2)  no  DNA  and  no  lipofectamine.  Both  controls  showed 
<0.4%  of  the  activity  measured  for  WT  HuGGT,  indicating  that 
the  GGT  activity  was  predominantly  that  of  the  recombinantly 
expressed  His-tagged  GGT  and  not  of  endogenous  GGT 
(a  membrane-associated  protein).  Protein  concentrations  of  the 
lysates  were  comparable  for  each  well  (data  not  shown).  The 
steady  state  kinetic  parameters  with  the  L-Cys  acceptor 
substrate  were  measured  using  the  WT  HuGGT  lysate  ( Kms  = 
4  ±  1  f.iM;  JCma  =  700  ±  300  // M)  and  were  comparable  to 
those  measured  using  purified  WT  HuGGT.  Mutation  of  Lys- 
562  to  asparagine  or  glutamine  led  to  a  4-  or  9-fold  reduction  in 
hydrolase  activity,  respectively.  The  loss  of  activity  is  consistent 
with  a  role  for  Lys-562  in  electrostatic  catalysis  and/ or  product 
release.  In  the  presence  of  the  glycyl-glycine  acceptor  substrate, 
a  57-  and  117-fold  decreases  in  transpeptidase  activity  were 
observed  for  the  K562N  and  K562Q,  variants,  respectively. 
The  decreases  in  transpeptidase  activity  that  were  larger  than 
those  in  hydrolase  activity  suggested  a  role  for  Lys-562  in 
neutralizing,  binding,  and/or  orienting  the  acceptor  substrate  in 
the  transpeptidation  reaction. 

Tyr-403  was  predicted  to  be  within  hydrogen  bonding 
distance  of  the  glycyl  group  of  the  glutathione  donor  substrate 
(Figure  6).  The  Y403A  variant  showed  approximately  3-  and 
14-fold  decreases  in  hydrolase  and  transpeptidase  activity  (with 
the  glycyl-glycine  acceptor  substrate),  respectively,  suggesting  a 


minor  role  in  catalysis  of  the  hydrolysis  reaction  and  a  more 
significant  role  in  the  transpeptidation  reaction. 

Binding  Interactions  of  the  Acceptor  Site  Inhibitor 
OU749.  To  further  probe  the  acceptor  site,  we  explored  the 
binding  interactions  of  OU749,  the  first  reported  acceptor  site 
inhibitor  of  a  GGT.34  OU749  is  a  potent  inhibitor  of  HuGGT; 
however,  it  does  not  inhibit  CapD  significantly  (IC50  >  1000 
I'M),  implying  that  residues  involved  in  binding  of  OU749  in 
the  acceptor  site  are  not  conserved  in  CapD.  Figure  7 A  shows 
the  effects  of  the  K562N,  K562Q,  and  Y403A  variants  on 
OU749  inhibition.  All  three  variants  significantly  increased  the 
IC50  from  30  ±  10  to  >2000  //M.  Consistent  with  the 
experimental  results,  the  predicted  binding  model  of  OU749 
with  HuGGT  showed  that  Lys-562  pointed  to  the  acceptor  site 
and  formed  a  hydrogen  bond  with  the  phenylcarboxylate  group 
of  OU749  (Figure  7B).  On  the  other  hand,  Tyr-403  formed 
aromatic  stacking  interactions  with  the  core  thiadiazole  group 
and  shielded  the  molecule  from  solvent.  Substitution  of  Lys- 
562  with  Gin  or  Asn  or  Tyr-403  with  Ala  disrupts  the  optimal 
binding  interaction  of  OU749  in  the  acceptor  site,  and  the 
calculated  binding  affinities  were  significantly  decreased. 

Interestingly,  the  phenyl  group  of  OU749  was  found  to  bind 
in  an  interior  pocket  of  the  putative  acceptor  site  of  HuGGT, 
which  was  not  observed  in  the  apo  structure  or  the  substrate 
binding  complex.  The  pocket  is  located  in  the  region  behind 
the  catalytic  residue  Thr-381  consisting  mainly  of  hydro- 
phobic  residues  such  as  Ala-471,  Ala-472,  Ala-542,  and  Val-544 
(Figure  7).  The  two  small  flexible  loops  [LI  and  L2  (Figure  5)] 
identified  from  MD  simulations  form  the  wall  at  the  bottom  of 
the  pocket.  In  comparison,  such  a  dynamic  hydrophobic  pocket 
was  not  observed  in  the  putative  acceptor  site  of  CapD.  Ser-427 
and  Pro-248  in  CapD  are  positioned  in  the  same  site  as  Ala-471 
and  Ala-472.  Docking  and  MD  simulations  showed  that  OU749 
did  not  fit  into  the  acceptor  site  of  CapD;  Ser-427  together 
with  Arg-520  appeared  to  hinder  OU749  binding. 

A  feature  of  OU749  inhibition  is  that  the  acceptor  site 
inhibitor  shows  high  selectivity  toward  HuGGT,  but  very  poor 
inhibition  of  rat,  mouse,  and  pig  GGT.34  The  species-specific 
inhibition  by  OU749  remains  an  interesting  issue,  as  these 
GGTs  are  highly  homologous,  and  the  structural  models  of 
HuGGT  and  rat  GGT  do  not  reveal  much  difference  in  the 
donor  or  acceptor  binding  sites.  Examination  of  the  OU749— 
HuGGT  complex  indicated  that  Ala-472  was  positioned  at  the 
entrance  of  the  hydrophobic  pocket  described  above  and  could 
play  a  role  in  regulating  access  to  this  binding  pocket.  Ala-472 
in  HuGGT  is  replaced  with  a  polar  Ser  in  rat,  mouse,  and  pig 
GGT.  Substitution  of  Ala-472  with  Ser  would  introduce  a  steric 
hindrance  into  the  hydrophobic  pocket.  Moreover,  the  side 
chain  could  form  a  hydrogen  bond  with  the  lysine  homologous 
to  Lys-562.  This  interaction  would  be  expected  to  block  the 
binding  of  OU749.  To  confirm  this,  we  modeled  one  HuGGT 
variant,  A472S,  and  performed  MD  simulations  to  probe  the 
dynamic  properties  of  the  hydrophobic  pocket.  As  shown  in 
Figure  7,  Ser-472  formed  a  stable  H-bond  with  Lys-562, 
resulting  in  a  closed  conformation  of  the  hydrophobic  binding 
pocket  in  the  acceptor  site,  whereas  in  WT  HuGGT,  Lys-562 
exhibited  high  flexibility  and  the  hydrophobic  pocket  remained 
in  an  open  conformation.  Therefore,  the  decreased  inhibitory 
activities  of  OU749  against  rat,  mouse,  and  pig  GGT  are  most 
likely  due  to  the  hydrogen-bonded  conformation  in  the  closed 
form,  which  would  prevent  OU749  binding. 
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Figure  7.  (A)  IC50  values  measured  with  WT  HuGGT  and  the  K562N,  K562Q,  and  Y403A  variants.  (B)  Binding  interactions  of  OU749  with 
HuGGT  (surface  representation).  The  OU749  inhibitor  and  the  active  site  residues  of  HuGGT  are  shown  as  sticks.  (C)  Open  conformation  of  the 
acceptor  site  of  HuGGT.  The  gating  residues  Lys-562  and  Ala-472  are  shown  in  a  surface  representation  with  the  carbons  colored  yellow,  oxygens 
red,  and  nitrogens  blue.  The  glutamate  substrate  (cyan)  and  catalytic  residue  Thr-381  (magenta)  are  shown  as  sticks.  (D)  Closed  conformation  of 
the  acceptor  site  observed  in  a  model  of  an  A472S  mutant.  The  gating  residue,  Lys-562,  can  form  an  H-bond  with  Ser-472.  (E)  Plot  of  distances 
between  residues  Lys-562  and  Ala-472  (Q,  atoms)  of  HuGGT  observed  in  the  MD  simulations.  (F)  Plot  of  distances  between  residues  Lys-562  and 
Ser-472  (C„  atoms)  of  the  A472S  mutant  during  the  MD  simulations. 


■  DISCUSSION 

Within  the  acceptor  site,  we  identified  a  positively  charged 
residue  (R520  in  CapD  and  K562  in  HuGGT)  that  is 
important  for  catalysis  and  acceptor  substrate  binding.  Despite 
differences  in  the  inhibitor  and  donor  substrate  preferences, 
both  CapD  and  rat  GGT  (79%  identical  to  HuGGT)  retain  the 
same  acceptor  substrate  preference  (both  prefer  L-Cys)  and 
have  the  lowest  Km3  for  L-Met  (Figure  2  and  Table  l).  While 
the  two  enzymes  can  utilize  the  same  acceptor  substrate,  they 
differ  in  their  susceptibilities  to  the  HuGGT  acceptor  site 
inhibitor,  OU749,  a  species-specific  inhibitor  that  does  not 
significantly  inhibit  CapD.  The  two  GGTs  also  differ  in  their 
acceptor  substrate  stereospecificity.  Rat  GGT  has  strict 


stereospecificity  toward  its  acceptor  substrates  and  can  utilize 
only  l- amino  acids  and  dipeptides;13’26  CapD  shows  no 
significant  preference  for  l-  or  D-amino  acid  acceptors  (Table  l) 
and  can  accommodate  both  stereoisomers.  Moreover,  CapD 
comparably  prefers  negatively  charged  acceptor  substrates  such 
as  d-  or  L-glutamate  when  compared  to  L-Cys.  From  the  site- 
directed  mutagenesis  and  modeling  described  herein,  it  is  clear 
that  Arg-432  and  Arg-520  play  an  important  role  in  acceptor 
binding  and  likely  allow  the  bacterial  enzyme  to  utilize  both 
d-  and  L-amino  acid  acceptors  via  two  alternative  binding 
modes  (Figure  3):  Arg-432  would  be  predicted  to  bind  the 
a-carboxylic  acid  of  the  amino  acid  for  a  D-isomer,  while  Arg-520 
would  bind  it  in  the  L-isomer.  EcGGT  has  not  been  shown  to 
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utilize  D-amino  acid  acceptors,27  and  this  may  be  due  to  its  lack 
of  an  arginine  equivalent  to  Arg-520  (EcGGT  has  a  serine  at 
this  position).  Correspondingly,  HuGGT  has  a  Gin  at  the 
position  equivalent  to  Arg-432  of  CapD  and  a  Lys  at  the 
position  equivalent  to  Arg-520  in  CapD  and  has  been  shown 
to  utilize  only  L-acceptor  substrates.  Several  bacterial  species 
( B .  anthracis,  F.  tularensis,  and  Y.  pestis )  have  GGTs  that  contain 
both  arginines  and  may  correspondingly  lack  acceptor  substrate 
stereospecificity. 

Mutation  of  Arg-520  led  to  the  most  significant  decrease  in 
fccat,  a  95-fold  decrease  for  hydrolysis  and  a  1400-fold  decrease 
for  transpeptidation.  The  modeling  and  kinetic  analysis 
suggested  a  role  for  Arg-520  as  an  electrostatic  catalyst.  The 
Arg-432  and  Arg-520  residues,  as  well  as  the  transiently 
protonated  N-terminal  amino  group  of  Thr-352  (which  has 
been  proposed  to  function  as  the  general  base  and  acid  in  Ntn- 
amidohydrolases14’21),  create  a  highly  positively  charged 
environment  in  the  acceptor  substrate  site.  Typical  acceptor 
substrates  would  place  a  zwitterion  near  Thr-352.  The  neutral 
amine,  however,  is  the  likely  nucleophile  in  the  deacylation 
reaction  (Scheme  SI  of  the  Supporting  Information).  The 
proximity  of  Arg-520  to  the  a-amino  group  of  the  acceptor 
substrate  would  likely  make  protonation  of  the  amino  group 
unfavorable  and  allow  the  neutral  amine  to  function  as  a 
nucleophile  in  the  transpeptidation  reaction.  Additionally,  the 
electrostatic  interaction  between  the  acceptor  COO-  group  and 
Arg-520  appears  to  be  important  for  acceptor  substrate  binding. 
Han  et  al.  speculated  that  a  positively  charged  residue  was 
present  within  the  acceptor  substrate  site.  On  the  basis  of  our 
kinetic  analysis,  modeling,  and  MD  simulations,  Lys-562  in 
HuGGT  and  Arg-520  in  CapD  appear  to  function  as 
electrostatic  catalysts  and  play  an  important  role  in  acceptor 
substrate  binding. 

The  transpeptidase  activity  has  been  proposed  to  serve  an 
important  physiological  function  and  accounts  for  at  least  50% 
of  the  utilization  of  the  y-glutamyl  donor.51  GGT  present  on 
the  surface  of  some  cell  membranes  is  important  for  the 
breakdown  of  extracellular  glutathione.  The  components  are 
utilized  or  recombined;  thus,  both  the  hydrolase  and  trans¬ 
peptidase  activity  appear  to  be  important  in  GSH  metabolism. 
In  B.  anthracis,  both  the  transpeptidase  and  hydrolase  activities 
function  in  anchoring  and  remodeling  the  capsule,  respectively. 
From  our  results  for  CapD  presented  here  and  the  previously 
published  results  for  rat  GGT  and  human  GGT,13'30  it  is 
apparent  that  the  transpeptidase  activity,  regardless  of  the 
amino  acid  acceptor,  is  higher  than  the  hydrolase  activity 
(Table  1  and  Figure  6).  Mechanistically,  it  has  been  unclear 
why  the  transpeptidase  activity  is  significantly  greater  than  the 
hydrolase  activity.  Examination  of  the  putative  acceptor 
substrate  site  suggests  that  the  negatively  charged  COO- 
group  of  the  acceptor  substrate  may  neutralize  the  positive 
charge  of  the  adjacent  arginines,  Arg-520  and  Arg-432,  and 
shield  the  a-amine  of  Thr-352  from  Arg-520.  Water  alone 
would  not  be  sufficient  to  neutralize  the  charge  of  Arg-520.  The 
shielding  likely  makes  the  protonation  of  the  a-amine  of  Thr- 
352  more  favorable  during  the  catalytic  cycle,  accounting  for 
the  faster  rate  of  transpeptidation  versus  hydrolysis  (Scheme  SI 
of  the  Supporting  Information). 

The  dynamics  of  the  acceptor  site  and  correlations  with  the 
donor  substrate  site  are  intriguing  characteristics  of  GGT. 
Experimental  studies  revealed  that  the  mechanism  of  inhibi¬ 
tion  of  HuGGT  by  OU749  was  uncompetitive  relative  to 
the  y-glutamyl  substrate  and  suggested  that  the  inhibitor 


preferentially  interacted  with  the  acylated  y-glutamyl  enzyme 
intermediate  over  the  free  enzyme.34  The  effect  of  the  y-glutamyl 
donor  on  the  acceptor  site  appears  to  be  subtle  as  the 
y-glutamyl  group  did  not  form  interactions  with  OU749  in  the 
HuGGT— OU749  complex  (Figure  7).  Analysis  of  the  dynamic 
behavior  and  conformational  changes  associated  with  the  apo 
and  y-glutamyl  enzyme  intermediate  in  HuGGT  indicated  that 
two  small  loops,  which  form  a  hydrophobic  pocket,  exhibited 
high  flexibility  in  the  substrate-free  form  but  were  significantly 
stabilized  in  the  glutamyl  enzyme  intermediate  complex  (Figure  5). 
The  hydrophobic  pocket  in  the  acceptor  site  of  HuGGT  was 
observed  in  the  substrate-bound  complex  but  not  apparent  in 
the  apo  state.  Therefore,  a  plausible  explanation  for  the 
uncompetitive  inhibition  of  OU749  is  that  the  formation  of  the 
covalently  bound  glutamyl  enzyme  intermediate  in  the  donor 
site  facilitates  rearrangement  of  the  acceptor  site  and 
consequently  favors  the  binding  of  acceptor  site  inhibitors  (or 
acceptor  substrates  during  transpeptidation). 

Our  overarching  goal  was  to  identify  distinct  differences 
between  the  bacterial  and  human  GGT  acceptor  sites  that  may 
be  useful  in  identifying  species-specific  inhibitors  and/ or  broad 
spectrum  antimicrobials  of  bacterial  GGT.  Here  we  demon¬ 
strated  that  the  use  of  small  molecule  probes  in  combination 
with  kinetic  and  computational  analyses  is  an  efficient  approach 
to  studying  the  underlying  structural  basis  of  substrate  and 
inhibitor  specificity.  The  characterization  of  the  acceptor 
substrate  site  of  HuGGT  and  B.  anthracis  CapD  provides 
insight  into  the  acceptor  and  donor  substrate  preferences 
associated  with  different  species.  More  importantly,  it  has 
revealed  a  dynamic  cooperation  between  the  donor  and 
acceptor  sites  in  substrate  recognition  and  catalysis.  This 
synergy  between  the  sites  may  be  useful  in  the  discovery  of 
potent  species-specific  bisite  inhibitors  that  bind  to  both  the 
donor  and  acceptor  sites. 
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or  /-glutamyl  transferase;  HpGGT,  H.  pylori  GGT;  HuGGT, 
human  GGT;  EcGGT,  E.  coli  GGT;  FtuGGT,  F.  tularensis 
GGT;  GSH,  glutathione;  MD,  molecular  dynamics;  PDGA, 
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